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We present a search for excited and exotic muon states /i*, conducted using an integrated lumi- 
nosity of 371 pb~^ of data collected in pp collisions at = 1.96 TeV at the Tevatron with the 
CDF II detector. We search for associated production of /x^* followed by the decay /i* fiy, 
resulting in the fijj,y final state. We compare the data to model predictions as a function of the mass 
of the excited muon M^* , the compositeness energy scale A, and the gauge coupling factor /. No 
signal above the standard model expectation is observed in the ^7 mass spectrum. In the contact 
interaction model, we exclude 107 < < 853 GeV/c^ for A = ; in the gauge-mediated 

model, we exclude 100 < M^. < 410 GeV/c^ for //A = 10"^ GeV"^ These 95% confidence level 
exclusions extend previous limits and are the first hadron collider results on jj,* production in the 
gauge-mediated model. 

PACS numbers: 12.60.Rc, 13.85.Qk, 12.60.-i, 14.60.Hi 



In the standard model (SM) the quarks and leptons of composite quarks and leptons consisting of fewer ele- 
are treated as fundamental particles. Their generational mentary particles than contained in the SM 0]. In this 
structure and observed mass hierarchy motivate a model model, quarks and leptons are the lowest-energy bound 
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states of these hypothetical particles, and additional ex- 
cited states exist near the compositeness energy scale A. 

Exotic fermions are also predicted in the context of 
grand unified or string theories, in which the known 
forces are unified into a larger symmetry group In 
such models, additional fermions are predicted with prop- 
erties similar to those of excited fermions. 

At a pp collider, excited or exotic muon states could be 
observed through the reaction qq (i* jjL. Excited muon 
production can be described using a contact interaction 
(CI) Lagrangian density [l|: 

47r _ 

L = -^qLl'^qLMLjf_,fiL + h.c, 

where Ml represents the left-handed fi* field. For exotic 
muon production, the relevant gauge-mediated (GM) La- 
grangian density is 0: 



L 



2A 



h.c 



where W and B are the SU{2)l and C/(l)y field 
strengths, g and g' are the respective electroweak cou- 
plings, and / and /' are phenomenological constants 
which are set equal to each other by convention to max- 
imize the photonic decay branching ratio. For / — 
f (/ = — /'), the relative /i* branching ratios are 
BR(^* ^ ^7) w 0.3 (0), BR(^* ^ vW) « 0.6 (0.6), and 
BR(^* ^ fiZ) w 0.1 (0.4) for M^, > 200 GeV/c^. The 
BR(^* fij) increases to 70% at M^. = 100 GeV/c^ 
for f — f- We use these branching ratios for both the 
GM and CI production models 

In this Letter we describe a search for associated fx^i* 
production that extends existing fj,* mass limits. Prior 
searches for /z* production have been performed by the 
LEP experiments, which have excluded fj,* with M^. < 
200 GeV/c^ for //A > 10"^ GeV'^ in the GM model 0. 
The D0 experiment has excluded n* with M^. < 688 
GeV/c^ for A = M^. in a particular CI model 0. 

We use 371 pb'^ of data collected with the CDF II de- 
tector at the Fermilab Tevatron. The detector has been 
described in detail elsewhere jS^ . Its magnetic spectrom- 
eter consists of silicon microstrip and drift chamber 
tracking detectors. Surrounding this are central @ (|?7| < 
1 0) and forward (hi > 1-1) electromagnetic (EM) 
and hadronic calorimeters. Embedded in the central EM 
calorimeter are wire and strip chambers (used to mea- 
sure the transverse shower profiles of photons) and a 
central preshower detector jy| (used for detecting pho- 
ton conversions). Outside the calorimeters are CMUP 
(|r?| < 0.6) and CMX (0.6 < \r/\ < 1) muon detectors jlll |. 
The momentum resolution of beam-constrained drift- 
chamber tracks is Spr/Px ~ 0.05%/(GeV/c). The elec- 
tromagnetic energy resolution for the high-energy pho- 
tons typical of fi* decays is approximately 2.5%. 

We analyze events passing the readout trigger require- 
ment of one drift-chamber track with px > 18 GeV/c |3 



matched to a reconstructed track segment in the muon 
chambers. In the offline analysis, we require two muon 
candidates identified by drift-chamber tracks with px > 
20 GeV/c and |?7| < 1, which pass requirements on 
impact parameter and number of hits, have minimum- 
ionizing particle properties, and at least one of which 
has a matching muon chamber segment [l^ . Both tracks 
must pass the calorimeter isolation requirement / < 0.1, 
where / is the ratio of the combined EM and hadronic 
calorimeter Et surrounding the track in a cone of radius 
R = ^(Ary)2 + (A(/))2 < 0.4 to the track pr- In ad- 
dition, both tracks must also pass the tracking isolation 
requirement T < 0.1, where T is the ratio of the summed 
Pt of tracks to the muon pT, where the sum is taken over 
tracks in a cone of radius i? < 0.4 around the muon di- 
rection. Finally, we reject cosmic rays based on tracking 
and track-timing information |l3l | . 

We select dimuon events that have a central or for- 
ward photon candidate with Ex > 25 GeV, \ri\ < 2.8, 
and / < 0.1, where calorimeter isolation is now defined 
relative to the photon. The photon energy is defined as 
the EM calorimeter energy in a 0.2 x 0.25 rj — cj) region, 
and its transverse component is determined by the vertex 
of the highest momentum muon with a matching muon 
chamber segment. Photon candidates are identified on 
the basis of their longitudinal and transverse calorime- 
ter shower profiles. Central photon candidates require 
the sum px of charged particles within a cone of radius 
R < 0.4 around the photon direction to be less than 
5 GeV + 0.005 x E'^, where E'^ is the Ex of the pho- 
ton, and no single charged particle within R < 0.4 to 
have PT > (1 + 0.005 x El/GeW) GeV/c. To suppress 
the initial-state radiation (ISR) Z{-^ background, 
we reject events with dimuon invariant mass m^^ in the 
range 81 - 101 GeV/c^. 

A Z — > /i/z sample is used to measure the efficien- 
cies of the muon identification criteria and trigger. The 
efficiency of the calorimeter and tracking identification 
requirements is measured to be (92.6 ± Q-2>stat)%- We 
measure the combined trigger and muon chamber match- 
ing efficiency to be (79.3 ± 1.05tat)% for the CMUP and 
(95.0 ± Q.&stat)% for the CMX. 

The photon identification efficiency is extracted from 
a GEANT-based detector simulation [lj|. Since the elec- 
tromagnetic showers of electrons are similar to those of 
photons, we validate the simulated photon efficiency us- 
ing a control sample of electrons from Z ^ ee events. 
The events are triggered and identified using one elec- 
tron; the second electron is used to emulate the pho- 
ton to determine the cut efficiencies for photons. In 
the central calorimeter, the emulated photon efficiency 
is (79.0 ± G.^stat)% in the data and (77.0 ± 0.&stat)% 
in the simulation. The corresponding values in the for- 
ward calorimeter are (89.0 ± 0.6stQt)% in the data and 
(89.6 ± Q-Qstat)% in the simulation. These comparisons 
show that photon detection is well-modeled by the simu- 
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lation; the statistically significant differences are included 
in the systematic uncertainty. 

The geometric acceptance and detector response are 
calculated with the GEANT detector simulation se par ately 
for the CI and GM models. We use the PYTHIA ^^ gen- 
erator for the CI model and the lanhep |3| and COM- 
PHEP ^3 programs to generate GM model events. The 
total signal acceptance (including identification efficien- 
cies) for the CI model increases from 13% at Af^. = 



to an asymptotic value of 21% for 



> 



100 GeV/c2 

400 GeV/c^ . The acceptance for the GM model increases 
from 12% at M^. = 100 GeV/c^ to an asymptotic value 
of 23% for M^. > 300 GeV/c^ . The relative systematic 
uncertainty on the acceptance is 3.1%, which is domi- 
nated by the uncertainty in the identification efficiency 
(2.2%) and simulation statistics (2.0%). 

We compute the expected background contributions 
from the following sources: (1) Z/'j* (2) 
Z/"f* (+7) — » Tr7, where the r's decay to muons, (3) 
Z/7* {—f fJ-i-J.)+ jet, where the jet is misidentified as a pho- 
ton, (4) t(— > iivb) -|-t(— > iJ-vb), where a fermion radiates a 
photon, (5) ev) -f Z{-^ fifi), where the electron is 

misidentified as a photon, and (6) Z(—f ee) -I- Z{-^ ^J■^J■)^ 
where one of the electrons is misidentified as a pho- 
ton. Other backgrounds (> 3 jets, W+ > 2 jets, and 
> 1 jet) were found to be negligible. 

The Zj, tt, WZ and ZZ backgrounds are estimated 
using simulated events, with the ZGAMMA generator 
for the Z{-^ l^-tAl background and pythia for the others. 
Systematic uncertainties on these background predictions 
arise due to integrated luminosity (6%) 0, parton dis- 
tribution functions (PDFs) (5%), higher-order QCD cor- 
rections (5%) [i^, acceptance (1%), and identification 
efficiencies (2%). 

The Z+jet background is estimated by weighting 
Z+jet events from the data by an i^T-dependent jet 7 
misidentification rate. The latter is measured using a 
jet-triggered data sample, correcting for the fraction of 
true prompt photons in the jet sample Q. The prompt 
photon fraction is estimated using 7 ^ ee conversions 
identified with the calorimeter preshower detector [2ll |. 
The jet^ 7 misidentification rate is applied as a func- 
tion of Et in the central calorimeter and as a function of 
Et and 77 in the forward calorimeter. 

We observe 17 signal candidates with a background 
prediction of 8.3 ± 0.9, of which 8.1 ± 0.8 are expected 
from Z"f production. The Poisson probability for the 
background to fiuctuate up to the observation, or higher, 
is 0.8%. Investigating further, we find that 11 candidate 
events have a 3-body mass m^;_j-y in the 81-101 GeV/c^ 
range, consistent with being final-state radiation (FSR) 
Z /i/i7 events, to be compared to the FSR Z /i/i7 
prediction of 5.5 ± 0.5 events. Figure □ and Table n show 
the m^j and m^^j distributions for the data and back- 
ground. 

To test our background prediction, we lower the mini- 



I I Z(— >|in,— >Tx)+Y 
r~1 Z+jet 
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FIG. 1: The fiy (a) and fifi-y (b) mass distributions for the 
data with background expectations. The total number of ob- 
served (expected) n^"/ entries is 17 (8.3 ±0.9). Both ^7 com- 
binations per event are included in (a). 



mum photon Et to 15 GeV and observe 43 events with a 
prediction of 38.5±4.0 events. The data show good agree- 
ment with expectation in this higher statistics sample, as 
shown in Figure |21 Additionally, we find consistency in 
the ISR Z7 control region of 81 < nifj,^ < 101 GeV/c^ 
and 25 < E^ < 50 GeV, where we observe 5 events with a 
prediction of 7.2^0 § events. We conclude that our signal 
sample has an upward statistical fluctuation, dominantly 
in the number of FSR Z fx^j events. 
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FIG. 2: 



The distributions for the data (43 events) and 
background (38.5±4.0 events) expectations for i5J > 15 GeV. 
The arrow indicates the minimum required for this anal- 
ysis. 



For the /j,* resonance search, we scan the m^^ spec- 
trum with a sliding window of width 3ct, where a is the 
mass width predicted by the simulation. Over almost the 
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TABLE I: Comparison of data and integrated background 
predictions above a given cut on the invariant mass of all ^'y 
combinations (left) and on the /1/17 invariant mass (right). 



a. 





combinations 




Events 




Data Background 




Data Background 


(GeV/c^ 


) 




(GeV/c2 ) 






> 


34 


16.6 ± 1.8 


> 


17 


8.3 ±0.9 


> 50 


22 


10.4 ± 1.1 


> 100 


6 


2.7 ±0.3 


> 100 


4 


2.1 ± 0.3 


> 150 


2 


1.5 ±0.2 


> 150 


2 


0.89 ±0.14 


> 200 


2 


0.9 ±0.1 


> 200 





0.37 ±0.07 


> 250 


1 


0.51 ±0.09 








> 300 





0.29 ±0.06 



entire model parameter space, a is dominated by detec- 
tor resolution. The tracker momentum scale and reso- 
lution, and the calorimeter energy scale and resolution, 
are tuned on the well-known Z — > /z^ and Z ee mass 
peaks j2^], respectively. For Af^» = A, the reconstructed 
/i7 mass resolution ranges from 9-90 GeV/c^ for masses 
ranging from 200-800 GeV/c^. 

We use a Bayesian ^22] approach, with a flat prior 
on the signal cross section and gamma priors on accep- 
tance and backgrounds, to set limits on the fifi* produc- 
tion cross section as a function of M^. The cross 
section limits are converted to M* limits by comparing 
them to the next-to-next-to- leading-order (NNLO) theo- 
retical cross sections j20] , computed using the MRST set 
of PDFs Hi]. We use the CTEQ prescription 26] to cal- 
culate the cross section uncertainty due to PDFs, which 
varies from 5% (M^. = 100 GeV/c^ ) to 20% (Af^- = 1 
TeV/c^ ). Uncertainties on higher-order QCD corrections 
(7-13%) depend on Af^- and the production model. 

Figure |3| shows the 95% confidence level (CL) upper 
limits on the experimental cross section x branching ratio 
and the theoretical curves for Af^. = A// (Af^. = A) 
for the CM (CI) model. Masses below 221 GeV/c^ (853 
GeV/c^) are excluded for the GM (CI) model for the case 
Af^. = A// (Af^. = A). Figure m shows the limits in 
the parameter space of / /A (A/^-» / A) versus Af^» for the 
GM (CI) model. For both models, we lose sensitivity for 
Mfj_* < 100 GeV/c^ due to large backgrounds and loss of 
signal acceptance in this range. We note that our CI limit 
of Mf,. > 853 GeV/c^ reduces to M^, > 696 GeV/c^ if 
we use /i* branching ratios that account for hypothetical 
CI decays, as assumed by the D0 collaboration 

In conclusion, we have presented a search for excited 
and exotic muons. In the GM model, we exclude 100 < 
Af^. < 410 GeV/c2 for //A = IQ-^ GeV"! at the 95% 
C.L., which is the first GM result at a hadron collider 
and is substantially beyond previous limits We have 
also presented ^* limits in the CI model as a function 
of Af^. and A, excluding 107 < Af^. < 853 GeV/c^ for 
A = A/m* . These limits complement our recent results of 
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FIG. 3: The experimental cross section x branching ratio 
limits at 95% CL for the CI (dashed-dotted hue) and GM 
models (solid line), compared to the CI model prediction for 
A — Mf_i* (dotted line) and the GM model prediction for 
A// = Af,t* (dashed line). Also indicated are the mass values 
that are excluded by these data. 
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FIG. 4: The 2-D parameter space regions excluded by this 
analysis for (a) the GM model, along with the current world 
limits 3, and (b) the CI model. We do not consider the 
region F^* > 2Af^* where the total width is larger than the 
11* mass. 
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